IFIC/12-81 
CFTP/12-019 



Exploring the Earth matter effect with 
atmospheric neutrinos in ice 



Sanjib Kumar Agarwalla a , TraceyLi b , Olga Mena^Jand Sergio Palomares-Ruizj^ 



(N 



O a,b,c i ns m u i de Fisica Corpuscular, CSIC-Universitat de Valencia, 

q Apartado de Correos 22085, E-46071 Valencia, Spain 

q d Centro de Fisica Tedrica de Particulas (CFTP), Instituto Superior Tecnico, 

Universidade Tecnica de Lisboa, Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal 



Abstract 

D 

We study the possibility to perform neutrino oscillation tomography and to determine the 
neutrino mass hierarchy in kilometer-scale ice Cerenkov detectors by means of the ^13-driven matter 
effects which occur during the propagation of atmospheric neutrinos deep through the Earth. 
OO We consider the ongoing IceCube/DeepCore neutrino observatory and future planned extensions, 
^vq such as the PINGU detector, which has a lower energy threshold. Our simulations include the 
impact of marginalization over the neutrino oscillation parameters and a fully correlated systematic 
uncertainty on the total number of events. For the current best-fit value of the mixing angle #13, 
the DeepCore detector, due to its relatively high-energy threshold, could only be sensitive to 
^ fluctuations on the normalization of the Earth's density of Ap ~ ±10% at ~ 1.6a CL after 
.£h 10 years in the case of a true normal hierarchy. For the two PINGU configurations we consider, 
^ overall density fluctuations of Ap ~ ±3% (±2%) could be measured at the 2a CL after 10 years, 
^ also in the case of a normal mass hierarchy. We also compare the prospects to determine the 
neutrino mass hierarchy in these three configurations and find that this could be achieved at the 
5<t CL, for both hierarchies, after 5 years in DeepCore and about 1 year in PINGU. This clearly 
shows the importance of lowering the energy threshold below 10 GeV so that detectors are fully 
sensitive to the resonant matter effects. 
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1 Introduction 



Cosmic ray interactions in the atmosphere are a natural source of neutrinos with baselines that 
span three orders of magnitude and energies in the range from MeV to well above TeV. This 
implies that different scales and effects could be accessible by studying atmospheric neutrinos. In 
fact, the most important result of the latest years in neutrino physics was the discovery in 1998 
of neutrino oscillations at the Super-Kamiokande detector by using atmospheric neutrino data in 
the GeV energy range [I]. 

Among the potential effects that could be explored are the resonant effects in the propagation 
of GeV atmospheric neutrinos through the Earth that mainly affect the subleading — > v e 
and v e — > Ufj, (or — > v e and v e — > v^) transitions. When atmospheric neutrinos deeply cross the 
Earth's mantle, the Mikheyev-Smirnov-Wolfenstein (MSW) resonance [2|3] could be in action j^-j6]. 
On the other hand, when they also traverse the core, resonant effects, different from the MSW 
resonance, could show up [7-15 . These effects are well known and their physics potential in iron- 
magnetized calorimeter, water-Cerenkov and liquid argon detectors has been thoroughly studied 
in the literature 



16 37 



It is very well known that the propagation of atmospheric neutrinos with energies in the 
range of a few GeV deeply through the Earth is sensitive to the currently unknown neutrino mass 
hierarchy [4| [l~6|[l~8|[2l] (or neutrino mass ordering), i.e., whether the third neutrino mass eigenstate 
is the lightest or the heaviest mass eigenstate. In the latter case, the hierarchy is called a normal 
hierarchy (NH) and in the former case, inverted hierarchy (IH). However, the ability to measure the 
neutrino mass hierarchy exploiting the resonant oscillation phenomena of atmospheric neutrinos 
inside the Earth crucially relies on the value of the mixing angle #13 . Recent measurements of this 



angle from the Daya Bay 38 , RENO 39 and Double Chooz 40 reactor experiments, in addition 
to the long-baseline T2K experiment 41 , seem to indicate that #13 ~ 9° 42 -44 . This large value 
of #13 opens up the possibility of using the atmospheric neutrino fluxes not only to determine 
the neutrino mass hierarchy, but also to study features of the Earth's matter density profile, by 
exploiting the #i3-driven matter effects. 

There are three different ways that have been proposed to infer some information about the 
internal structure of the Earth by exploiting the different effects of neutrino propagation in matter: 
neutrino absorption tomography, neutrino oscillation tomography and neutrino diffraction. The 
idea of using the absorption of very high energy neutrinos to explore the Earth's interior dates 



back to 1974 45 and is analogous to X-ray tomography, but instead exploits the attenuation 
of the neutrino flux for energies E v > 10 TeV 46 . There are numerous studies which have 



considered neutrinos of different origins, such as man-made neutrinos 45 , 47 56 , extraterrestrial 



neutrino^] [5~T| [57|^60] and atmospheric neutrinos 61 64 . On the other hand, neutrino oscillation 
tomography relies on the matter effects in neutrino oscillations and it has been considered by 
studying man-made beams 65-75, solar 76,77 and supernova neutrinos 77,78. The third 



possibility is based on studying the diffraction pattern of crystalline matter in the interior of the 



1 This possibility was first suggested by J. Learned and H. Bradner in the late 1970s, but we lack a proper 
reference. 
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Earth caused by coherent neutrino scattering, but it is technologically unfeasible 79 . 



In the case of atmospheric neutrinos, only neutrino absorption tomography has been considered 

However, resonant effects are very sensitive 



within the context of kilometer-scale detectors 61 - 64 



to the matter density that neutrinos traverse, so one could also think of doing neutrino oscillation 
tomography with atmospheric neutrinos with energies in the GeV rang fl This is one of the main 
goals of this work and represents a completely different technique to determine the Earth's density 
distribution to those used in geophysics, mainly based on the analysis of seismic waves. Although, 
in principle, geophysics can obtain more precise results, its inferences are based on numerous 
assumptions because the velocities of seismic waves not only depend on the density, but also on the 
composition, temperature, pressure and elastic properties of the medium. Thus, the results from 
atmospheric neutrino tomography would represent an independent and complementary assessment 
of the Earth's internal structure. 

Neutrino telescopes such as AMANDA 



80 



lected a large number of atmospheric neutrino events 83|-|89 
threshold 



IceCube [81] and Antares 82 have already col- 
even in spite of their high energy 
100 GeV) and the steeply falling atmospheric neutrino spectrum (~ E~ 3 ' 7 ). Whereas 



this is a very high threshold for studying neutrino oscillations with atmospheric neutrinos, a new 
generation of neutrino telescopes with lower energy thresholds (~ 10 GeV), such as the DeepCore 
extension of the Icecube detector, is currently taking data successfully [90}]92] and further natural 
extensions of this are being planned. The Precision IceCube Next Generation Upgrade (PINGU) 
has been proposed in order to reduce the detection threshold down to a few GeV (93] . Although 
reaching those energy thresholds in these multi-Mton scale neutrino telescopes is very challenging, 
if successful, the atmospheric neutrino events detected at them would offer a great opportunity for 
detailed oscillation studies, including the determination of the neutrino mass hierarchy 94,95 , tan 



neutrino appearance searches 96 , and precise measurements of the atmospheric neutrino oscilla- 



tion parameters 97 . Given the large value of (9 13 measured by the different reactor 38 -40 and 



the T2K 41 experiments, the prospects look promising. It therefore seems timely to assess their 



ability to infer some information about the Earth's matter density profile and how this could also 
affect the determination of the neutrino mass hierarchy from the atmospheric neutrino events. In 
this work we study the future prospects for DeepCore and PINGU by considering the experimental 
capabilities of the two detectors and present the expected sensitivities to both of these observables. 

The structure of the paper is as follows. We start in Sec. [2] by revisiting the atmospheric 
neutrino fluxes, the Earth's internal structure and the matter effects which affect the propagation 
of GeV neutrinos inside the Earth. Sec. [3] contains a description of the detectors considered here, 
as well as of the x 2 analysis methods used to calculate the sensitivities to the Earth's matter 
density and neutrino mass hierarchy, which are presented in Sec. |4j Finally, in Sec. [5] we draw our 
conclusions. 



2 This was suggested in Ref. 
distributions were studied. 
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, where the changes in the transition probabilities for different Earth density 
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2 Atmospheric neutrino fluxes and oscillations in the Earth 



2.1 Atmospheric neutrino flux 



When cosmic rays traverse the Earth's atmosphere, hadronic showers are produced by their inter- 
actions with the particles that constitute it. Depending on the energy of the primary cosmic ray, 
different secondaries can be produced. Below ~ 100 GeV, the neutrino flux is dominated by the 
pion decay chain, whereas above these energies, kaon decays dominate neutrino production. The 
decay chain is the following: 



p + X -»■ 



Ti- 



ll 



--/K± 



Y 



— > ji z 



(1) 



// ± — > e* + v e 



v,. 



Thus, when the conditions for the decay of all these particles are satisfied, we expect flavor ratios 
for the atmospheric neutrino fluxes of v^fv^ ~ 1 and (z/ M + Vy)j{v e + z/ e ) ~ 2. This flux is usually 
referred to as the conventional atmospheric neutrino flux and the energies of these atmospheric 
neutrinos range from a few MeV up to hundreds of TeV. 



The gross features of the spectrum are easy to understand (see e.g., Refs. 98 991). For energies 



E v <C 100 GeV/| cos#|, with 6 being the zenith angle, all pions decay and the neutrino spectrum 



has approximately the same power as the primary cosmic-ray spectrum, which is ~ E~ . Above 
these energies, the expected power-law neutrino spectrum (from pion decays) is asymptotically 
one power steeper. This has to do with the extra factor of 1/E n in the ratio of the pion decay 
length to interaction length. The fact that the flux is higher in the horizontal direction (it is 
proportional to sec#) is explained by the longer decay path near the horizon: a pion traveling 
through the atmosphere horizontally (cos# ~ 0) has a higher probability of decaying in the 
atmosphere than a pion traversing the atmosphere vertically (cos 6 ~ 1) and therefore the flux is 
higher for the horizontal component. A similar behavior occurs for kaons, but shifted to energies 
about one order of magnitude higher. On the other hand, the muon spectrum is very similar to 
the parent meson spectrum, but for energies above ~ 1 GeV, many muons reach the ground before 
decaying (at ~ 100 GeV, only ~ 15% of atmospheric muon neutrinos come from muon decay), and 
hence the neutrino flavor ratio grows with increasing energy. Similarly to the angular dependence 
of neutrinos from pion and kaon decays, this results in the steepening of the electron neutrino 
spectrum occurring at lower energies for vertical directions and at higher energies for propagation 
close to the horizon. Existing analytical approximations to the atmospheric neutrino fluxes versus 
energy and zenith angle at high energies also allow a good qualitative understanding of all these 



features 98,99 



The huge range of energies and baselines provided by the atmospheric neutrinos opens up 
the possibility of exploring many different and exciting physics topics, from the measurement of 
neutrino mixing parameters 100 101 to more exotic ones 84 102 114 . At very high energies, 



101 to more exotic ones 84 102 11 
above 10 TeV, neutrino interaction cross sections become high enough for neutrinos going through 
the Earth to start becoming attenuated. This effect is sensitive to the neutrino interaction cross 
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sections and to the density profile of the Earth and thus, a measurement of the atmospheric 
neutrino-induced event rate at these energies could provide information about the Earth's internal 
structure 61-64 . The intermediate-energy region, between ~ 50 GeV and ~ 1 TeV, provides 
information about the atmospheric neutrino flux normalization 84 . Finally, in the low-energy 



region which we exploit here (below about 50 GeV), neutrino oscillation effects could be sig- 
nificant [4,(7-16,18,21 , in particular due to the large value of the #13 mixing angle recently 
measured 38 -41 . Resonant matter effects in the propagation of neutrinos inside the Earth could 
be very important for E v ~ [1, 15] GeV and could strongly enhance or decrease the oscillation 
probabilities, for neutrinos in the case of normal neutrino mass hierarchy and for antineutrinos if 
the neutrino mass hierarchy is inverted. 

The atmospheric neutrino fluxes we use in this work are those from Ref. 



115] (see also 



Refs. 116 118 for other atmospheric neutrino flux calculation^]). The absolute electron and 
muon atmospheric (anti)neutrino fluxes are uncertain at the level of 10% - 20% in the energy re- 
gion of interest here, which is related to our ignorance regarding hadron production models. The 
uncertainties quoted above are largely canceled for the muon neutrino-antineutrino flavor ratio as 
well as for the up-down ratio, the former expected to be of order 1% above 1 GeV 119 . Accel- 



erator data from the HARP, MIPP and NA61 experiments on particle multi-production are also 
expected to improve the predictions of the absolute atmospheric neutrino fluxes 120 . However, 



as we will see when discussing the effects of correlated systematic uncertainties in the analysis, 
the uncertainties related to the normalization of the atmospheric fluxes are alleviated when higher 
energy bins, as well as different angular bins, in which oscillation effects are negligible, are also 
used. 



2.2 The Earth's internal structure 

Most of the knowledge we have about the internal structure of the Earth and the physical properties 
of its different layers comes from geophysics and, in particular, from the data we obtain from seismic 
waves. Other information, from geomagnetic and geodynamical data, solid state theory and high 
temperature-pressure experimental results is also used. It turns out that a reliable estimate of the 
density of the Earth is essential to solve a number of important problems in geophysics, such as 
the dynamics of the core and mantle, the gravity field, the mechanism of the geomagnetic dynamo, 
the bulk composition of the Earth, etc. (see, e.g., Ref. |121| and references therein). 

The Earth is conventionally divided into three main (approximately) spherical concentric shells: 
the crust, the mantle and the core, each of them further divided into subshells, with different 
properties (see, e.g., Refs. 122 123] ). Whereas the crust only represents about 0.4% of the Earth's 



mass, the mantle and the core constitute about 68% and 32%, respectively. On the other hand, 
the core has a radius almost half that of the Earth and it is about twice as dense as the mantle. 



3 Note that in Ref. 118 the atmospheric neutrino flux at the South Pole was calculated. At this site there is no 
rigidity cutoff, and the expected flux is slightly higher than that at the Kamioka site where there is an intermediate 
rigidity cutoff. However, the main differences show up at low energies {E v <C 10 GeV) and are not expected to 
significantly affect our results. 
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The divisions of the Earth's interior are based on the reflection and refraction of compressional 
and shear waves, i.e., P- and S-waves, respectively. The crust is the brittle outer shell, with 
a large proportion of incompatible elements. It is composed of three different crustal types: 
continental (thick and composed primarily of granite), transitional (defined by the islands, island 
arcs and continental margins) and oceanic (thinner and composed primarily of basalt). The crust 
is separated from the mantle by a boundary discovered by Mohorovicic in 1909. The mantle 
is composed mainly of silicates and oxides, being a poor conductor of electricity and heat and 
with very high temperature-dependent viscosity. Its temperature increases with depth creating 
a geothermal gradient which is responsible for the different rock behaviors that determine the 
mantle subdivisions. Whereas in the upper mantle rocks are cool and brittle, in the lower mantle 
they are hot and soft, but not molten. In 1891 E. Wiechert inferred the existence of a completely 
different region beneath the mantle from the mean density and moment of inertia of the Earth, but 
it was not until 1914 that B. Gutenberg made the first accurate determination of the location of 
the core-mantle boundary 124 from the R. D. Oldham measurement in 1906 of a sharp decrease 
in the velocities of P- and S-waves deep in the Earth. Soon after this, the first models of the 
radial structure of elastic velocities started to be developed. This boundary presents a density 
contrast of a factor of ~ 2 and a viscosity contrast of ~ 20-24 orders of magnitude, which is what 
ultimately causes the strongly differing time and length scales for convection motions in these two 
layers. The Earth's core, the Earth's source of internal heat, is thought to be composed mainly of 
an alloy of iron and nickel, with a small admixture of lighter elements, which is a good conductor 
of electricity and heat. Its radius is approximately 2900 km and it is divided into the outer (liquid) 
and the inner (solid) cores. 

Although the Earth is not exactly spherically symmetric and is irregular, the use of mean mod- 
els, taking the Earth's physical parameters as symmetric, is very useful as a reference framework 



for different studies 125 . During many years, the most widely used radially symmetric model 
has been the Preliminary Reference Earth Model (PREM) |126 , which included anisotropy in the 
upper mantle, and it was a good fit to free oscillation center frequency measurements, surface- 
wave dispersion observations, travel-time data for body-waves and some astronomical data such 
as the Earth's radius, mass and moment of inertia. However, there was no discussion about the 
sources of error or covariances. A more recent 1-D model, AK135, which combines travel times 



and free oscillations, is also available 127,128 . Although there are observed deviations at the few 
per cent level that cannot be explained by 1-D models, to first order, the structure of the Earth 
is determined by its radial properties and these reference models are used as a starting point for 
more refined studies. In Fig. [I] we show the matter density profile of the Earth (and fluctuations 
of ±10%) as a function of the distance to the center according to the PREM |126| . Up to a 
radius of ~ 2900 km, the density increases gradually from the crust to the mantle. Then, in the 
core-mantle boundary, there is a very sharp transition, followed by a gradually increasing density 
until the center of the Earth. 

For the last half a century the problem of determining the density distribution of the Earth has 
mainly consisted of perturbation inversion using different seismic data. This is a very demanding 
problem, which in many cases is non-linear. Moreover, most studies of the Earth's radial structure 
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Figure 1: The Earth's density profile according to the PREM \12Q{ (red) and with a ±10% overall density 
fluctuation (green dash-dotted and blue dashed lines), as a function of the distance from the center of the Earth. 



make inferences about different properties given a specific density model or based on empirical re- 
lations between seismic waves velocities and density. Nevertheless, there are significant trade-offs 
with other crucial parameters because wave velocities also depend on composition, temperature, 
pressure and elastic properties, which necessarily imply some modeling and uncertainties in the 
density estimate. Recently, 3-D models, using 1-D models as a reference, have been developed, al- 
though some studies concluded that few robust density features could be constrained with existing 
data [129- 132] (see, however, Refs. 133-138]). All in all, density values, averaged over ~ 100 km 
are thought to be known at the level of a few per cent at all depths, whereas density gradients are 
less known. On the other hand, linear integral constraints on the spherically symmetric density 
distribution are known at the level of ~ 10~ 4 from estimates of the Earth's mass and moment of 
inertia 139 141] (see also Ref. 138 for a recent analysis of the spectra of the lowest frequency 
seismic mode with a precision of ~ 10~ 3 ). This implies that global variations of the density are 
constrained to be within those uncertainties. Obviously then, overall density variations at the 
level of a few per cent (as those shown in Fig. [l]) are already excluded. Nevertheless, we will not 
impose these linear integral constraints and evaluate, using the PREM as our reference model, 
how sensitive atmospheric neutrino oscillation tomography is to changes in the density. 



2.3 Oscillation probabilities 

For neutrino energies in the range of a few GeV, the transition probabilities v„ — > v e 



and v P v, 



— > v^t)) of atmospheric neutrinos in their propagation through the Earth are 
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relevant if genuine 3-flavor neutrino mixing takes place, i.e., if the #13 mixing angle is different from 
zero [4 15 . Moreover, in this energy range and for these baselines (L > 1000 km), CP- violation 



effects are very small and can be safely neglected. Likewise, effects due to the 1-2 sector are also 
subdominant and, as a first approximation, can also be neglected. In this context, the calculation of 
the transition probabilities effectively reduces to a 2-neutrino problem, with Am^ and #13 playing 
the role of the relevant 2-neutrino oscillation parameters. Within these approximations, the 3- 
neutrino oscillation probabilities of interest for atmospheric v e ^ having energy E v and crossing 
the Earth along a trajectory characterized by a zenith angle^j^, have the following form [7 11 13 
(see also Refs. §[lO|[T4||T5)): 

PtoiVe^Ve) - 1-P2u, (2) 
P$v{Ve -> Vy) ~ P*v(Vp U e ) ~ Sin 2 6*23 P2v , (3) 
Pto(v e ->V T ) ~ COS 2 9 23 Piu , (4) 

P 3 u{^ ^ 1 - \ sin 2 29 23 - sin 4 23 P 2v + i sin 2 29 23 Re {e~ iK A 2v {v T -+ v T )) , (5) 

P-iv^y. ~+V T ) = 1 - P 3u (Uy -> - P 3l/ (l/ A4 ->■ I/ e ) (6) 



where P 2v = P2 U (Aml l ,9i 3 ; E u ,9) is the 2-neutrino probability describing v e — > v x transitions, 
where v x = sin 823 + cos #23 v T , and k and A2 V [y T — >• i/ T ) = ^21/ are the phase and the 2-neutrino 
transition probability amplitude. For antineutrinos the oscillation probabilities are analogous to 
those for neutrinos: they can be obtained formally from Eqs. Q - (|6p by changing the sign of the 
matter potential (or equivalently, p by —p). It is interesting to note that, within the approximation 
Auijj = 0, the probabilities for neutrinos and NH (IH) are the same as those for antineutrinos 
and IH (NH). 

Therefore, the magnitude of the matter effects depends on the 2-neutrino oscillation probability 
Piv- In case of oscillations in vacuum, P 2v ~ sin 2 2#i3, so this probability is small. However, matter 
effects can strongly enhance P2 V and thereby greatly modify the 3-neutrino probabilities. On the 
other hand, if #13 = 0, then P 2u = and Re (e~ lK 'A 2v {y T — V u T )) — cos (Am^L/ (2 E v )), and hence 
matter effects are absent. If this were the situation, these probabilities would get reduced to the 
case of -H- v T 2-neutrino oscillations, so for NH and IH they would be equal and identical to the 
case of vacuum oscillations. 

In Fig.[2]we show the — > (upper panels) and v e — > (lower panels) transition probabilities 
as a function of cos 9 for two different energies (E v = 5 GeV in left panels and E v = 10 GeV in right 
panels) for NH (red regions) and IH (green regions limited by blue lines) . The bands correspond to 
how the probabilities change if the density (according to the PREM) varies up to ±10%. We can 
see that matter effects, that are very sensitive to the value of the density, tend to greatly enhance 
the v e -r-7- Up transitions and reduce the — > survival probability, and also shift the positions 
of the maxima and minima with respect to the case of negligible matter effects. Moreover, a 
change in the matter density profile would shift the location of both the resonance energy and 

4 Neutrinos with cos# = — 1 are directly upgoing and traverse the entire diameter of the Earth, those with 
cos 9 = come from the horizon and those with cos 9 > are downgoing and reach the detector from above. 
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COS0 COS0 

Figure 2: Oscillation probabilities. Upper panels: Piv^ as a function of cos 9, for E v = 5 GeV (7eft 

panel) and E u = 10 GeV (right panel), for NH (red regions) and IH (green regions limited by blue lines). Lower 
panels: Same but for P(v e — > v^). The widths of the bands correspond to varying the matter density by ±10%. 
Similar results are obtained for the case of antineutrinos by exchanging the curves for NH (IH) by those for IH 
(NH). We have used the current best-fit values for the oscillation parameters f42| . 
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the baseline at which the resonant behavior is maximized, modifying the number of events in 
the different angular and energy bins. For propagation in the mantle only, the resonance is, to a 
good approximation, the MSW resonance. However, in the case of propagation through the core 
(cos 9 < —0.83), non-trivial resonant effects show up, although for the energies considered in Fig. [2] 
they are not maximal. Notice that matter effects are larger for E v = 5 GeV and more baselines 
are affected than for E v = 10 GeV, the reason being that resonances are closer to E v = 5 GeV 
(see below). 

As we discussed above, the Earth is composed, in addition to the crust, of two major density 
structures that constitute almost all the mass of the planet: the mantle and the core. Although 
in our calculations we compute the full 3-neutrino evolution and use the detailed PREM for the 
density distribution^} as a first approximation and in order to gain insight about the physics 
of neutrino propagation in the Earth, one can consider a simplified model for the density profile, 
divided into two shells with different densities, p m ~ 4.5 g/cm and p c ~ 11.5 g/cm , with the core- 
mantle boundary at a radius of ~ 2900 km. All the interesting features of the atmospheric neutrino 
oscillations in the Earth can be understood quite accurately within this framework [7[jl5]. There 
are analytical solutions for the transition probabilities for neutrinos crossing the Earth [7f{l6|[l8 



but they reduce to the case of neutrino propagation in a medium of constant density for trajectories 
such that cos 9 > —0.83, i.e., for neutrinos which propagate only in the mantle and could experience 
MSW resonant effects In this case, Eqs. ((2]) — (JsJ) reac 



'Am j 

i <=;r, 2 oa m <=;r, 2 ' 31 
l — sin ^#13 sin 



'Am t 

sin 2 02a sin 2 29™ sin 2 ' 31 



'Am r 

cos 2 9^ sin 2 29™ sin 2 ' 31 



1 
2 



1 - r sin 2 2023 - sin 4 23 sin 2 29™ sin 2 ' 31 



sin 2 29 



23 



A 3 i L 
cos I — - — 1 — 



'13 

.4 



Am 2 ! 



cos 29™, sin 



A 3 iL 



.4 



Am 2 x 



sm 



A™L 



+ 



cos 



A™L 



A™L 



(7) 



(9) 
(10) 



5 Moreover, throughout this work we will consider as a good approximation the Earth to be an isoscalar neutral 
medium, i.e., n e /(n p + n n ) = 1/2, where n e ,n pi n n are the electron, proton and neutron densities, respectively. 

6 Note that for a spherically symmetric Earth, the neutrino trajectory is fully specified by its zenith angle, 
L = -2 Ra cos 9. 
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where 



Am?, , s 

A 3 i = , 12 

3 2E ' 



j i> 



AS = — yJ(Am 2 31 cos20 13 -A)2 + (Am| 1 sm2M 2 ) (13) 
sin 2 2^3 ee sin 2 20 13 (A > ) , (14) 



^A 

and the matter parameter is A = 2 y/2Gp n e E u , with Gp the Fermi constant and n e the electron 
number density. The expressions for antineutrinos can be obtained by replacing A — > —A. Thus, 
the resonant behavior, when maximal mixing occurs, happens for the case of NH (IH) in the 
neutrino (antineutrino) channel at the resonant energy, 

Am 2 cos 26 13 _ /4.5g/cm 3 \ / Am 2 \ nn M . 

E res ee ~ 7GeV - =■ cos 26>i 3 . (15) 

2V2G F n e \ P y V2.4 x 10- 3 eV 2 ; 1 ; 

On the other hand, the first oscillation maximum of the oscillation term happens when the oscil- 
lation phase is = |. Hence, the baseline at which both the conditions for the resonance and 
the first oscillation maximum are satisfied, is [2] 

W = * - 1-1 ' ^ ( 4M ^) ( ] ■ (16) 



\[2G F n e tan 26 13 \ P I V tan2 #i3 



For sin 2 26 Vi ~ 0.1, L max ~ 1.1-10 4 km and E Tes ~ 7 GeV for densities approximately corresponding 
to the average value of the mantle. However, strong matter effects are also present around the 
resonant energy even if the baseline does not correspond exactly to the first oscillation maximum, 
i.e., even for baselines L ~ L max /\/2 For illustration purposes, we show in Fig. [3] the line 
average density for a given trajectory in the Earth (upper panel) and the corresponding resonance 
energy in the case of propagation in a medium of constant density with that value (lower panel). 

Using Eqs. ^ - ^ it is straightforward to see that the fluxes of atmospheric v e>fi of energy E v 
which reach the detector after crossing the Earth along a given trajectory, ^ Ue ^(E u ,6), are given 
by the following expressions in the case of the 3-neutrino oscillations under discussion |9- Hj 18] : 

0) * < [l + (sin 2 fl 23 r-l)P 2 ,] , (17) 
- \ sin 2 2# 23 + 



^(E u ,6) ~ ^ 

sin 4 # 23 ((sin 2 #23 r)- 1 - l) P 2v + \ sin 2 2 6 23 Re {e~ iK A 2v {v T -»■ v T )) 
where $° = $° (E v , 6) is the v e (jj) H ux & t the production point in the atmosphere and 



(18) 



rs ^ s *&d>- (19) 
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Figure 3: Line average of the Earth's density and resonance energy. Upper panel: Line average 
density of the Earth for a given trajectory with total baseline L. Lower panel: Resonance energy for propagation 
of neutrinos in the Earth assuming the density is constant and equal to the line average density. We have assumed 
sin 2 = 0.025, Am^ = 2.44 x 10~ 3 eV 2 and the density distribution according to the PREM (red lines) and with 
a ±10% overall density fluctuation (green dash-dotted and blue dashed lines). 



For neutrino trajectories with cos# < —0.4, i.e., L > 5100 km, the predicted ratio of to 
v e atmospheric neutrino fluxes is r ~ 2.0 — 2.2 (2.1 — 2.6) for sub-GeV (anti) neutrinos while 
2.1 — 5.6 (2.4 — 7.2) for multi-GeV (anti)neutrinos [115 1 18] . Thus, even in the case of resonant 



r 



matter effects taking place, the effects of #i 3 -driven transitions are suppressed for the sub-GeV 
sample due to the term (sin 2 #23 r — 1), especially if sin 2 #23 ^ 0.5 (the current 3<r confidence level 



(CL) allowed range is sin 2 623 = (0.36 — 0.68) 42 ). On the other hand, for multi-GeV neutrinos, 
(sin 2 023r — 1) > 0, so actually this factor can amplify the matter effects. Nevertheless, neutrino 
telescopes are Cerenkov detectors, so they have no charge-identification capabilities and therefore 
cannot distinguish neutrinos from antineutrinos. As the resonances take place for NH for neutrinos 
and IH for antineutrinos, the ^13-driven matter effects would consequently be smeared and the 
sensitivity to them reduced as compared to the case when measurements of the neutrino-induced 
rates and antineutrino-induced rates can be performed separately However, notice that the 
fact that neutrino and antineutrino cross sections are different could allow us to distinguish NH 

7 The distinction between the neutrino and antineutrino channels can be achieved with magnetized detectors 
such as the magnetized iron calorimeter at the India-based Neutrino Observatory (INO) |142 143 , but with much 
smaller volumes than neutrino telescopes. 
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from IH. Unfortunately, neutrino telescopes have very poor angular resolution for cascades, so the 
electron neutrino-induced event rates, where these resonant effects are expected to be larger (see 
Fig. [2j, would give little information. Therefore, muon tracks are the main signal in this study 
and hence, the muon neutrino-induced event rates. 

Following the discussion above, from Eqs. ^ - Eqs. (17) - (18) and the analogous equations 
for antineutrinos, we see that ^-driven matter effects increase the <H- v e (and z/ M <H- u e ) 
transitions with sin 2 #23, leading to an increase of the electron neutrino-induced event rate and 
a decrease of the muon neutrino-induced event rate at the detector. The rates are larger for 
multi-GeV events than for sub-GeV events. 



3 Analysis 

3.1 Detectors set-up 



The IceCube/DeepCore detector 90 is a densely instrumented region located at the bottom 
center of the IceCube detector at a depth of between 2100 m and 2450 m, in such a way that 
it avoids a horizontal layer with a high content of dust and therefore poor optical properties. 
The detector consists of six additional strings instrumented with phototubes of higher quantum 
efficiency with respect to IceCube. Thus, in general, DeepCore has multiple advantages when 
compared to Icecube: at the depth at which DeepCore is located, ice is on average twice as clear 
as the ice above, allowing for a lower neutrino energy threshold. Additionally, the larger amount 
of photosensors, separated by 7 m instead of 17 m for IceCube, and the higher quantum efficiency, 
lead to a significant gain in sensitivity of up to a factor of 6, especially for low energy neutrinos. 
The rest of the IceCube detector, along with a horizontal region with additional instrumentation 
at a depth of 1800 m, could be used as an active veto for downgoing atmospheric muons, allowing 
the study also of downgoing atmospheric neutrinos. 

At the energies of interest here, the angular reconstruction capabilities of DeepCore for showers 
is not optimal. Therefore, for the present study, we only consider upgoing muon-like events with 
the effective volume for the 86-string configuration (IC86) at trigger (SMT3) and online filter level 
in the 10-15 GeV energy region as shown in Tab. [I] |144| . We follow a very conservative approach 
by only considering this single energy bin (E v = [10, 15] GeV) and discarding lower energy bins 
which in principle have a lower effective volume and worse detection capabilities. The IceCube 
Collaboration aims to achieve a signal efficiency of ~ 50% for contained and partially contained 



events 90 



There is also a recent proposal to further upgrade the IceCube detector, with the planned 



PINGU 93 . It would consist of 20 additional strings within the DeepCore volume to lower the 
neutrino energy threshold down to 0(1) GeV energies and the detector would also benefit from 
the DeepCore strings. 

In this work, we consider two possible PINGU scenarios: a conservative scenario which we refer 
to as PINGU-0 with 5 GeV energy threshold and with two 5 GeV energy bins [E v = [5, 10] GeV 
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Detector 


Energy bins (GeV) 


p ice V eS (Mton) 




so in 


3 91 


l J L. L, | ) V LJ I 


1 fl 875 - 1 9 85 


4 46 




1 2 8,5 - 14 895 


5 69 




14 825 - 16 8 


6 53 




^ 6 


4 91 




6 7 


4 87 




7-8 


4 Q1 




8 - 9 


5.79 


PINGU 


9 - 10 


7.54 




10 - 11 


6.73 




11 - 12 


6.77 




12 - 13 


7.58 




13 - 14 


7.86 




14 - 15 


8.77 



Table 1: Detectors set-up. Effective mass for DeepCore and PINGU for the relevant energy bins which are 
considered in the current study (see text for details) \144\ . In all cases, we add a post-trigger detection efficiency 
of 50% (not included in these numbers). For DeepCore we consider a single energy bin, E v = [10, 15] GeV, and 
for PINGU we consider events in the interval E v = [5, 15] GeV but we use two different sets of energy bins: for 
PINGU-0, E v = [5, 10] GeV and E v = [10, 15] GeV; and for PINGU-I, E v = [5, 7.5] GeV, E v = [7.5, 10] GeV, 
E v = [10, 12.5] GeV and E v = [12.5, 15] GeV. The size of the angular bins is always 0.1 in cos9 and we consider 
events in the interval cos 9 = [—1, 0]. 

and [10, 15] GeV), and a less conservative scenario which we call PINGU-I, with four 2.5 GeV 
energy bins [E v = [5.0, 7.5] GeV, [7.5, 10.0] GeV, [10.0, 12.5] GeV and [12.5, 15.0] GeV). The 
effective masses for these configurations (Tab. [I]) are computed assuming a trigger setting of 
3 digital optical modules hit in 2.5 ps that are in local coincidence, a containment criterium 
which is implemented by a cut on the ^-position that matches the DeepCore fiducial volume 
(—500 m < z < —157 m) and a tight radius cut from string 36 (r < 150 m) which is the center of 
DeepCore/PINGU [lii] . 

Throughout this work we assume that the determination of the neutrino energy and direction 
could be achieved with some degree of accuracy. This implies that in addition to the characteri- 
zation of the muon tracks induced by the atmospheric muon neutrino interactions in the detector, 
the accompanying hadronic cascades must also be detected. This is a challenging task, so in this 
respect, the results presented here should be taken as optimistic. Nevertheless, we have assumed 
very wide energy bins, so that a neutrino energy resolution at the level of < 20% is not expected 
to significantly affect our results. On the other hand, if the angular resolution is not better than 
the typical root mean square value of the scattering angle between the incoming neutrino and the 
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produced muon, V ^ ~ yjGeV / E u , some smearing is expected and hence, the sensitivity would 
be slightly worse. Let us also note that the effective masses used in this work represent preliminary 
MonteCarlo results based on the outcome of ongoing simulations by the IceCube Collaboration 
and these numbers might be further refined in the near future. 

In summary, for the DeepCore detector simulations we assume a single energy bin of E u = 
[10, 15] GeV. For the PINGU-0 scenario, we assume two energy bins: E u = [5, 10] GeV and 
[10, 15] GeV. For PINGU-I, we assume four energy bins: E v = [5.0, 7.5] GeV, [7.5, 10.0] GeV, 
[10.0, 12.5] GeV and [12.5, 15.0] GeV. The effective masses we consider are shown in Tab. [T] In 
all configurations we assume that the post-trigger efficiency of the detector, for all energy bins, is 
50% and consider an angular binning in cos# of width 0.1 for cos# G [—1, 0]. 

3.2 Number of muon-like events 

Cerenkov detectors do not have charge-identification capabilities, so neutrinos cannot be distin- 
guished from antineutrinos. Hence, these two types of events must be added together in the 
analysis. The number of contained atmospheric muon-like events (neutrino- plus antineutrino- 
induced samples) at the detector in a given angular (i) and energy bin (j), is given by 

r(cos6) i+1 ME u ) j+1 

(N^j = 2nN A T / dcosO dE v e p ice V eS (E„) x 

J{cos8)i J{Ev)j 



d cos 9 dE v d cos 6 dE v 



\dcosudEy dcosvdEp J 



(20) 



where is the Avogadro number, T is the time of data taking, e = 0.5 is the post-trigger 
efficiency, p lce is the density of ice, V e s is the effective volume (the effective mass, p- lce V e g, is given 
in Tab. [I] for the different configurations we consider) and <j v qq and are the deep inelastic 
scattering cross sections for neutrinos and antineutrinos off an isoscalar target^ 

a v c N c ~ 7.30 ■ 10- 39 (E„/GeV) cm 2 ; (21) 
a D c N c ~ 3.77 ■ 10- 39 (E„/GeV) cm 2 . (22) 

The number of atmospheric neutrino- and antineutrino-induced contained events that are ex- 
pected in the PINGU-0 set-up after an exposure of 10 years is shown in Fig. [4] as a function of 
cos#. The left panel corresponds to the 5-10 GeV bin and the right panel to the 10-15 GeV 
energy bin. The number of events for both neutrino mass hierarchies, NH and IH, are shown for 
neutrinos and antineutrinos. We also depict two cases where we increase the Earth's density by 



8 Note that ice is not an isoscalar target. However, the differences in the cross sections with respect to an isoscalar 
target are at the level of 2-3%. Notice also that the values we use are slightly different from the world-averaged 



values over energies that extend well above 100 GeV 145 ■ 147 . 
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Figure 4: Number of muon-like contained events in PINGU-0 after 10 years for NH and IH in the 

case of a PREM matter profile, in the 5-10 GeV bin (left) and 10-15 GeV (right) as a function of cos 9 (also shown 
is the corresponding baseline on the upper axis). We also show the event numbers in the case of a PREM matter 
profile with an overall fluctuation of +10% for neutrino with NH and antineutrino with IH. Note that the vertical 
axis is different in each panel. Although shown independently, in the analyses we add together the neutrino and 
antineutrino events, due to the lack of charge-identification of the detectors. 



10% with respect to the PREM: neutrino-induced events for NH and antineutrino-induced events 
for IH, where resonant matter effects are at play. We show the neutrino and antineutrino rates 
independently for the sake of illustrating the matter effect taking place for one case or the other, 
depending on the mass hierarchy, although in the analysis they are added. 

These plots show the dependence of the numbers of muon-like events on the mass hierarchy 
(compare the red solid and the blue dotted lines for neutrinos with the magenta dotted and the 
black dotted lines for antineutrinos) and on changes of the normalization of the matter profile 
(compare the red solid with the green dotted lines for neutrinos and NH and the magenta dotted 
with the brown dotted lines for antineutrinos and IH). We only show the results with a change 
in the density with respect to the PREM for the neutrino-NH events and the antineutrino-IH 
events because matter effects take place inside the Earth in the case of NH for antineutrinos and 
of IH for antineutrinos. The effect is larger for the neutrino channel mainly due to the larger cross 
section and due to the fact that the flux of neutrinos is slightly higher than that of antineutrinos, 



Vfxl ~ 1.15 and v e /v e ~ 1.25, for the energies of interest here 115-118 



A general feature in both panels is the absence of significant matter effects for cos 6 > —0.4 (L < 
5100 km), as expected (see Sec. 2.3). For neutrinos crossing the Earth deeply, resonant 6^3-driven 
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matter effects are very important and there are clear differences in the number of events between 
neutrino-NH and neutrino-IH (antineutrino-NH and antineutrino-IH) , and between neutrino-NH 
(antineutrino-IH) for a density profile given by the PREM and for the case of increasing the 
density by 10%. In particular, notice that there are important differences in the total number of 
muon-like events in the angular interval —0.8 < cosd < —0.6 (7645 km < L < 10194 km) for 
the two cases of the depicted Earth's density, being always more important for NH as it is in this 
case when the neutrino channel is mostly affected by matter effects. The larger the density the 
lower the resonant energy and the shorter the baseline of the first oscillation maximum. Hence, 
the crossing of the number of events of these two cases (e.g., the red solid and green dotted lines) 
for —0.6 < cos8 < —0.5 is due to the interplay between these two effects. On the other hand, 
the matter effect on the propagation of core-crossing neutrinos is smaller than for just-mantle- 
crossing neutrinos. This is due to the resonance energy for the former being outside the energy 
range we consider here, and hence the effect is due to non-resonant matter effects that modify the 
propagation in a less significant manner. 

The matter effect is larger in the lower part of the energy range we consider, as is clear from 
both panels and was discussed above. In general, the resonance energy for just-mantle-crossing 
neutrinos lies within the 5-10 GeV energy bin, whereas in the 10-15 GeV energy bin the effects 
are non-resonant for all trajectories. The minimum at cos 6 ~ —0.5 for 10-15 GeV (right panel) 
is explained by the minimum in the P{y v — > z/ M ) probability (upper-right plot in Fig. [2]), which 
moves to smaller values of cos 9 (longer baselines) for higher energies. Likewise, the minimum 
at cos# ~ —0.3 for 5-10 GeV is explained by the first minimum in the P{y v — > Vy) probability 
(upper- left and upper-right plots in Fig. [2]). The fact that the second minimum in this probability 
does not translate into a minimum in the number of muon-like events in the energy range 5-10 GeV 
is because the first maximum contributes significantly as the full energy range is integrated. 

Overall, the best sensitivity to the neutrino mass hierarchy (difference between the red solid 
and blue dotted lines, and between the magenta dotted and black dotted lines) and to the matter 
density (difference between the red solid and green dotted lines, and between magenta dotted and 
brown dotted lines) is achieved in the 5-10 GeV energy range and for neutrino trajectories such 
that —0.8 < cos# < —0.4, although some sensitivity remains for other angular and energy bins. 
In this case, the resonant matter effects take place, whereas for higher energies or for neutrinos 
traversing the core, non-resonant effects occur. As expected, it is clear from both figures that the 
sensitivity to the mass hierarchy is better than that to changes in the Earth's density. 

3.3 Statistical analysis 

We have generated our simulated data for the true (central) values of the oscillation parameters 
given in the first column of Tab. [2] The numbers of simulated muon-like events in each bin for 
10 years of data taking in DeepCore and PINGU-0 are given in Tabs. [4] and [5] in Appendix |A} We 
also show in these tables the results after varying the normalization of the matter density. 

The values of the oscillation parameters we use to estimate the sensitivities lie within the 
allowed ranges at la CL obtained from global fits of the current neutrino data, except for the 
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True Values 


Marginalization Range 


External la error 


sin 2 ^ ue = 0.025 


[0.019, 0.0301 


a(sin 2 29 n ) = 5% 


sin 2 #* r o ue = 0.5 


[0.38, 0.661 

L i J 


cr(sin 2 2 9 23 ) = 2% 


(Am 2 ff ) tmc = ±2.4 x lO" 3 eV 2 


[2.2, 2.6] x 10~ 3 eV 2 (NH) 
-[2.6, 2.2] x 10~ 3 eV 2 (IH) 


a(|Am 2 ff |)=4% 


(Am 2 21 ) truc = 7.62 x 10" 5 eV 2 






sin 2 9\ 2 nc = 0.32 






<%p e = 0° 






Ap tmc = 


[-0.1, 0.1] 




ttruc n 

Snorm u 


[-1, 1] 


^norm 20% 



Table 2: Benchmark parameters used in the analyses. The first column shows the central values of the 
oscillation parameters used in this work as true values. The second column shows the ranges of the parameters 
over which the marginalizations have been performed in the Ht (2a CL ranges '421). Estimated future 1<t CL 
errors on the oscillation parameters are given in the last column. In the last two rows, Ap is the correction to the 
normalization of the Earth 's matter density with respect to the PREM prohle and £ norm represents the nuisance 
parameter used to incorporate a fully correlated systematic error (a norrn ) in the normalization of the number of 
events. 



value of sin 9 



23 



42 



44 



Recent global analyses prefer a non-maximal value of 9 2 % at about the 

although the new 



2a CL which is mainly due to the recent MINOS z/„ disappearance data 148 



Super-Kamiokande atmospheric neutrino data still favors maximal 2-3 mixing as the best-fit value 

0.5 as the 



in a 2-neutrino oscillation analysis 149 



In our analysis, we choose to use sin 2 9\ rm 



'23 



reference value. For the atmospheric neutrino mass splitting, we take an effective mass-squared 
difference of (Am 2 ff ) tme = ±2.4 x 10 -3 eV 2 as measured by the accelerator experiments in the 

The + (— ) sign refers to the case where NH (IH) is the true 



v, 



— > Vn disappearance channel 148 



hierarchy. This effective mass-squared difference is related to the Am^ and Am^ mass-squared 



differences through the expression 150 151 



Am 2 ! 



Am 21 (cos — cos<5r:p sin^?. sin26 l 1 5 tan# 



'12 



CP 



'13 



'12 



y 23j 



(23) 



where Am^ = mf — rn 2 -. Since the degenerate solutions occur at ±Am 2 ff , which is slightly different 
from ±Am 2 ;1 , it is important to take into account this feature when evaluating the sensitivity to 
the neutrino mass hierarchy. 

In all fits, we marginalize over sin 2 26 l i3, sin 2 2#23, and |Am 2 ff | within their presently allowed 
± 2a ranges (see the second column of Tab. [2]). We impose a prior (Gaussian error at la) of 5% on 
sin 2 2#i3 considering the expected precision by the Daya Bay experiment by 2016 152 . We also 
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take expected priors of 2% and 4% on sin 2 29 2 3 and |Am 2 ff |, respectively |153| . We keep 6*12, Am 
and Sqp fixed both in data and in theory and no marginalization has been performed for these 
oscillation parameters since they would affect our results in a negligible way. Finally, we also add 
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a fully correlated systematic error in the normalization of the number of events, which could be 
due to errors on the normalization of the atmospheric neutrino flux, the detector effective mass, 
the cross section or the efficiency We set this error to a noTm = 0.2. 

The treatment of correlated systematic errors is performed by the Lagrange multiplier method 

Thus, we consider nuisance parameters that describe the systematic 



or pull approach 154-157 



and the errors of sin 26 l i 3 (£, 



sin 2 2#* r 3 uc ) 



and |Am 2 ff | (£, 

Am , 



sin 2 20i3 

|Am 2 ff | 



error of the normalization of the number of events (£„ 
sin 2 2^3 - sin 2 2^), sin 2 2fl 23 (& m 2 26 , 23 = sin 2 2# 23 

\(Am 2 s ) trnc \). The variation of these nuisance parameters in the fit is constrained by adding a 
quadratic penalty to the corresponding \ 2 function without systematic errors. 

Let the number of muon-like events detected in the i-th angular and j-th energy bin be 



Ardata 
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(24) 



with h = sign(Am 2 ff ) and Ap a global normalization factor for the Earth's density so that p = 
(1 + Ap) Pprem; with Pprem being the PREM density profile. We estimate the difference between 
data and theory with the following \ 2 expression 154 157 : 



X 2 (h, Ap) =min (5normiCsin2 , 2 , 



Xstat norm I 



+ 



Snorm 



+ X 



prior 



(25) 
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and Xprior IS defined as 
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(27) 



In order to estimate the sensitivity to the global normalization of the Earth's matter density, 
in the simulated experimental data we set Ap tme = 0, and in the fit we vary Ap in the range 
[—0.1, 0.1]. Then, we marginalize over the two possibilities for the neutrino mass hierarchy (h = 
±), so that we obtain x 2 (Ap). The sensitivity is calculated as 

S(Ap) = y/ X 2 (Ap) - x 2 (0) = y/A X 2 (Ap) . (28) 

We present our results for both NH (h tme = +) and IH (h tine = — ) as the true mass hierarchy. 
In the case of the determination of the neutrino mass hierarchy, we define 



xl(Ap) = min 



a 2 201 



e 23 '€[Am 2 ff |) 



XmH±(A=F' £norm) + 



(Tu 



+ X 



prior 



(29) 
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for true NH (x^(Ap)) and true IH (x_(Ap)), where 

(cose), iV ij ; 

with A + (A_) for NH (IH) defined as 

A± = {013, 023, |Am 2 ff |, ±, Ap; 6\™, (Am^) ta , 5*™} . (31) 

Then we marginalize over Ap in the range [—0.1, 0.1] to incorporate the impact of possible fluctu- 
ations in the global normalization of the density of all the layers of the PREM profile (but without 
imposing any prior) to obtain xl- The sensitivity in this case is just given by 

Sk = ^ h . (32) 

Like for the sensitivity to the Earth's density, we also study both cases for the neutrino mass 
hierarchy, when data is generated assuming NH (IH) as the true neutrino mass hierarchy and we 
compare with the theoretical expectation when considering IH (NH). 

Notice that for both analyses, only the few bins where the #i3-driven matter effects are present 
(see Fig. [i] and Tabs. |4] and [H] in Appendix [A]) would contribute to the \ 2 - Thus, correlated 
systematic errors are expected to have a much smaller effect than uncorrelated systematic errors 
(see Appendix [B]), as the values of the parameters related to them get fixed by the bins where the 
matter effect is negligible. We have explicitly checked this for both analyses. 



4 Results 

In this section we show the sensitivities to the Earth's density distribution and to the neutrino 
mass hierarchy by using atmospheric neutrino data for the DeepCore, PINGU-0 and PINGU-I 
configurations. We also comment on the impact of systematic errors on the analyses, as well as 
the effect of the marginalization over the oscillation parameters |Am 2 ff |, #13 and #23. 

4.1 Sensitivity to the Earth's matter density 

In Fig. [5] we show our results for the sensitivity to the Earth's matter density, for both NH 
(left panel) and IH (right panel) as the true hierarchy. We show the sensitivities for the three 
detector configurations detailed in Tab. [T] for an exposure of 10 years. All the results shown 



in this section include the effects of the marginalization, within the current 2a CL ranges 42 
over sin 2 20x3, sin 2 26*23 an d |Am 2 ff |, on which we impose a prior as discussed in the previous 
section. We also marginalize over the neutrino mass hierarchy. Finally, we also include a 20% 
fully correlated systematic uncertainty in the total number of muon-like events. As mentioned 
above, this systematic error does not have a large impact on the sensitivities due mainly to the 



19 



25 



15 



< 



- 1 1 \: 

: NHtrue 


• • : 


: 

• • • 

■ • 


; • - 
■ • 
• • 






■ 

»• ■ 

**, 


, » » 




/ • 








* • 


**, 
*• 


PINGU-I J; ,* : 

/:/ - 


- .»• 

3a : *• 
■ ** 


<• * ; 

/ »*PINGU-0 : 


»• 
*• 

r 

~ ~2a~ " ? " ^ 


/ •*: 

»■« - - 

• - - /V- - - 

*• DfcepCore^^^ 



10 



< 



-10 



-5 





Ap [%] 



10 




-10 



. 1 1 1 1 1 1 1 1 1 
'- i IH true 


,,,,,,, 

« 




«f 


i 


*: 
• - 


\ : 3a 
• 


» - 
» - 

-, -i — 


: » 
: • 
» 

> \ 
• • 


P.INGU-I; : 
• 
• 
• 


^ « 


• 

5 
• 


* *, 


• * 


* * 

: — _»_. V 

i \\\ 2a 

V 

L \\ 

DeepCore ■ 

r* i i i Sto 


/ / 

// 
/ / 

/ ;*PINGU-0: 

iMV i i i — 



-5 





Ap [%] 



10 



Figure 5: Sensitivity to fluctuations on the overall normalization of the Earth's density with respect 
to the PREM profile, for the DeepCore (red solid lines), PfNGU-0 (blue dashed lines) and PfNGU-f (brown dash- 
dotted lines) conhgurations for the case of NH as the true hierarchy (left panel) and for fH as the true hierarchy 
(right panel). 

presence of the higher cos 9 (closer to the horizon and thus very little affected by matter effects) 
angular bins in the analysis, that help to fix the normalization of the number of events in those 
bins. 

For the DeepCore configuration, fluctuations in the normalization of the Earth's density of 
Ap ~ ±10% can be detected at the ~ 1.6a CL (~ lo~ CL) in the case where NH (IH) is the true 
hierarchy. The prospects are significantly improved for the PINGU-0 and PINGU-I configurations, 
due to the information contained in the lower energy bins (matter effects are maximal at the 
resonant energies ~ 5 — 7 GeV, which are unaccessible to the considered DeepCore configuration). 
For NH (IH), PINGU-0 would be able to measure matter fluctuations of Ap ~ ±3% (±9%) at the 
2d CL. PINGU-I would improve these numbers to Ap ~ ±2% (±6%) for NH (IH). 

The results are always significantly better for NH. This can be understood from the fact that 
the matter effects are present in the neutrino channel, which has more statistics. For IH the larger 
number of events in the neutrino channel (not affected by matter effects in this case) substantially 
swamp the ^13-driven matter effects in the antineutrino channel. As we have mentioned, these 
have to be added together, which reduces the sensitivity to density fluctuations, specially in the 
case of IH. On the other hand, it is interesting to note that PINGU-0 and PINGU-I are slightly 
more sensitive to negative fluctuations of the density than to positive ones. This is because the 
higher the density, the lower the resonance energy (see lower panel of Fig. [31). Thus, with an energy 
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threshold of 5 GeV, some of the events affected by the matter effects would lie outside the energy- 
range we consider. However, for lower densities, more events around the resonance region are 
included in the data set, which compensates the loss of sensitivity due the marginalizations. This 
has a smaller effect in DeepCore since the resonant effects do not take place in the energy bin we 
consider. We have also explicitly checked that, for the assumed energy thresholds, these detectors 
are basically sensitive to fluctuations of the normalization of the Earth's mantle (L < 10700 km), 
as the resonance for cross- crossing neutrinos occurs at lower energies. 

Finally, let us comment on the impact of the marginalization over all the different parameters, 
which is small. Indeed, for negative density fluctuations, it is negligible for PINGU-0 and PINGU- 
I. For positive fluctuations, the marginalization over 623 has a small but non-negligible effect for 
the three detector configurations. The value of sin 2 #23 sets the amplitude of the matter effects (see 
Eqs. ([2]) - (|6| and Eqs. ^ - (10)) and therefore, in principle, one expects its marginalization to 
have some impact. As a consequence of what was just discussed above, the effect would be larger 
for the case of positive fluctuations for which the sensitivity is slightly worse. The configuration 
for which the largest effects of the marginalizations occurs is PINGU-0 (at the level of < 10%). 
DeepCore does not include the resonance region and hence, the effects of the marginalizations are 
similar for negative and positive fluctuations. On the other hand, PINGU-I has more energy bins, 
which helps to distinguish density effects from changes in the mixing parameters. 

Notice that the measurement of the matter density profile is much more difficult than the 
extraction of the neutrino mass hierarchy, which we consider next. This can be explained by 
the fact that while just measuring matter effects alone can serve to determine the neutrino mass 
spectrum, for the Earth's matter density one also has to be able to measure perturbations to the 
matter effects. These could be small and difficult to disentangle unless the energy and angular 
information of the signal is reasonably good. 



4.2 Sensitivity to the neutrino mass hierarchy 

The sensitivity to the neutrino mass hierarchy by exploiting the number of muon-like events in the 
three possible neutrino configurations considered in this work is shown in Tab. |3j which contains 
the results of the \ 2 analyses for different exposure times (1, 5 and 10 years). We show the results 
for the two possible choices of the true neutrino mass hierarchy, NH and IH. As in the case of 
the study of the matter density fluctuations, for the analysis of the determination of the neutrino 
mass hierarchy, we also marginalize over the current 2a allowed range of the oscillation parameters 
sin 2 2#i3, sin 2 2^23 and | Am 2 ff |, imposing the same priors expected from future measurements, and 
add a 20% fully correlated systematic error related to the overall normalization of the number 
of events. In addition, here we also marginalize over the Earth's matter density distribution, by 
allowing its normalization to vary freely within ±10%. 

From the results in Tab. [3] we can see that the prospects of measuring the neutrino mass 
hierarchy with these three detector configurations are very promising. The ongoing DeepCore 
detector with the signal efficiency and the energy reconstruction capabilities assumed here could 
provide a 3cr (5a) measurement of the neutrino mass hierarchy after slightly more than 1 year 
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Detector configuration 


Exposure [years] 


CL [a] NH(true) 


CL [a] IH(true) 




1 


2.90 


2.63 






fi 17 

U. 1 1 


5 36 




10 


8.38 


6.92 




1 


4.98 


4.91 


PINGU-0 


5 


10.2 


8.49 




10 


14.0 


11.0 




1 


6.50 


6.22 


PINGU-I 


5 


13.3 


10.8 




10 


18.4 


14.5 



Table 3: Sensitivity to the neutrino mass hierarchy for the three detector configurations explored here for 
different exposure times. We show the results for NH (third column ) and IH (fourth column ) as the true hierarchy. 



(less than 5 years) if nature has chosen NH and after less than 2 years (5 years) for IH. As we 
have already mentioned, the results are worse for the IH scenario because in this case the resonant 
behavior occurs in the antineutrino channel, which is statistically suppressed due to the smaller 
antineutrino cross sections and initial fluxes. On the other hand, with detector configurations such 
as those considered here for PINGU-0 and PINGU-I, the mass hierarchy could be measured with 
an astonishing precision, ~ 5<r and > 6a after 1 year for both NH and IH, respectively. 



We have checked that our results agree with previous findings in the literature 94,95 . However 



note that those analyses differ on a number of important points from the one presented here. In 
particular, in the recent Ref. [95] , a larger effective volume for PINGU was used and no post-trigger 
efficiency was considered. Overall this amounts to about a factor of four more statistics in their 
case. In that work, a much finer binning, both in E v and cos#, was also considered. This is mainly 
the reason why they get a sensitivity to the mass hierarchy for 5 years of 45.5<r CL for NH for the 
case of no systematic errors and no smearing, whereas for PINGU-0 (PINGU-I), we get 10.2a CL 
(13.3a CL). On the other hand, whereas we consider large bins for the neutrino energy and zenith 
angle to approximately take into account the uncertainties in their reconstruction which smear 



the signal, in Ref. 95 these uncertainties are studied more explicitly. Even with the much larger 
bins we consider, we would expect our results to slightly worsen when adding this information in a 
more accurate way. Another important difference is the way in which systematic errors are added. 
Whereas they add uncorrelated errors on the normalization of the number of events, we add fully 
correlated uncertainties. As briefly discussed in Appendix [B] and checked numerically, this has 
important consequences when only a few out of all the considered bins contribute to the x 2 , as is 
the case here. In the case of uncorrelated systematic errors in the normalization, the sensitivity 
gets reduced as the systematic uncertainty increases, whereas for correlated systematic errors the 
sensitivity is not much affected with respect to the case of no systematic errors, thanks to the use 
of bins where this error gets fixed. Hence, their results get worse the larger the systematic error 
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is, whereas the impact of the systematic error on our results is very small. Finally, we also use 
different central values for the oscillation parameters, although this has a small effect. 

Finally, note that the neutrino mass hierarchy sensitivities computed here are almost indepen- 
dent of the value of the CP violating phase, #cp- Therefore the mass hierarchy determination for 
these three detector configurations would be a clean measurement, free of degeneracies, and could 



serve as an input for other experiments which focus on CP violation searches, such as the T2K 41 



or NOz/A |158| long-baseline neutrino experiments (see Ref. 159 for a recent analysis combining 
the capabilities of both experiments). 



5 Conclusions 

Atmospheric neutrinos produced by the interactions of cosmic rays in the Earth's atmosphere 
cover a huge range of energies and baselines, offering the opportunity to address many open 
physics questions. In this work we have exploited the low energy region, E u = 5 — 15 GeV, where 
matter effects could affect the propagation of these atmospheric neutrinos while they pass deep 
through the Earth. If # 13 is sufficiently large, the transition probabilities — > u e {p^ — > u e ) and 
^ e — > v ^(t) {pe — > ^/i(r)) of atmospheric neutrinos become relevant and resonant matter effects 
would affect the neutrino (antineutrino) oscillation channels if nature has chosen NH (IH). Hence, 



the recent measurement of #13 ~ 9° from several reactor neutrino oscillation experiments 38 -40 



as well as from the accelerator based T2K experiment 41 , brings the opportunity to study these 
#i3-driven matter effects. This could allow us not only to determine the neutrino mass hierarchy, 
but also to infer some general features of the Earth's density profile. 

Indeed, the idea of exploiting matter effects in atmospheric neutrino oscillations to distinguish 
the type of neutrino mass hierarchy has been extensively explored in the literature (4 37 , 94 



95 . On the other hand, three different ways to study the density profile of the Earth by using 



neutrinos have been considered: neutrino absorption tomography 45 , 47 - 64 , neutrino oscillation 



tomography 65 - 72 , 74 78 and neutrino diffraction 79 . These techniques have been discussed 
for accelerator, extra-terrestrial, atmospheric, solar and supernova neutrinos. With regards to 
exploiting the atmospheric neutrino flux, only neutrino absorption tomography had been studied in 
some detail. Here we have considered neutrino oscillation tomography with atmospheric neutrinos. 

Obtaining a reliable estimate of the density of the Earth is essential to solve a number of 
important problems in geophysics. In this work we have presented a completely different technique 
to determine the Earth's density profile to those used in geophysics, which are mainly based on 
the analysis of the velocities of seismic waves and on empirical relations between these velocities 
and density. Although, in principle, geophysics can obtain more precise results, its inferences 
are based on numerous assumptions as there are significant trade-offs with other parameters. 
Thus, the results from atmospheric neutrino tomography would represent an independent and 
complementary assessment of the Earth's internal structure. 

The goal of this work is the study of the Earth matter effects taking place for atmospheric 
neutrinos in the GeV range. We have explored the possibility to do neutrino oscillation tomography 
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to infer information on the Earth's matter density and have studied future sensitivities to determine 
the neutrino mass hierarchy. In order to do so, we have considered kilometer-scale ice Cerenkov 



detectors such as the ongoing DeepCore extension of IceCube 90 92 , and PINGU |93|, which is 



a further proposed extension to IceCube. Although these multi-Mton detectors have no charge- 
identification capabilities that would allow us to distinguish neutrinos from antineutrinos and 
hence increase the sensitivity to matter effects, the different cross sections for neutrinos and 
antineutrinos (and also the slightly different fluxes), the large statistics which can be accumulated 
by these detectors, and the large value measured for the mixing angle # 13 , make the study of 
these matter effects feasible and make possible the determination of the neutrino mass hierarchy 
in relatively short time scales. 

In this work, three different experimental scenarios have been examined (see Tab. |l|). For the 
simulations of the DeepCore set-up we have assumed a single energy bin of E u = [10, 15] GeV 
and an effective mass of ~5 Mton. For the proposed PINGU experiment, we have considered 
two possible configurations. For the PINGU-0 set-up, we have assumed two energy bins: E u = 
[5, 10] GeV and [10, 15] GeV and an effective mass of ~ 7 Mton. For the PINGU-I set-up 
we have assumed four energy bins: E v = [5.0, 7.5] GeV, [7.5, 10.0] GeV, [10.0, 12.5] GeV and 
[12.5, 15.0] GeV and an effective mass of ~ 7 Mton. In addition, in all simulations we have 
considered a signal efficiency of 50%. 

In our analyses we have addressed the impact of marginalizations over the neutrino oscillation 
parameters and we have also included fully correlated systematic uncertainties (20%) on the overall 
normalization of the total number of events. Our results are shown in Fig. [5] and Tab. |3j For 
the current best fit values of the oscillation parameters (see Tab. |2|, the DeepCore detector, with 
its relatively high-energy threshold, would only be sensitive to fluctuations on the normalization 
of the Earth's density of Ap ~ 10% at ~ 1.6a CL (la CL) after 10 years for NH (IH). On the 
other hand, the measurement of the neutrino mass hierarchy, being a discrete parameter, is much 
easier and DeepCore could provide a > 5a CL measurement after 5 years of data taking, for 
both hierarchies. For the PINGU detector, the lower energy thresholds assumed for the PINGU-0 
and PINGU-I set-ups provide better sensitivities to both the Earth's density fluctuations and the 
neutrino mass hierarchy. For NH, overall density fluctuations of Ap ~ ±3% (±2%) with respect 
to the PREM profile could be measured at 2a CL after 10 years in PINGU-0 (PINGU-I). In the IH 
case, the former numbers translate into Ap ~ ±9% (±6%). A measurement of the neutrino mass 
hierarchy at the > 5a CL might be possible within 1 year of data taking, for both hierarchies. This 
clearly shows the importance of lowering the energy threshold below 10 GeV, so that detectors 
are fully sensitive to the resonant #i 3 -driven matter effects. 

In summary, we have studied the future prospects and sensitivities for the ongoing DeepCore 
and the proposed PINGU detectors to infer overall fluctuations in the Earth's matter density 
profile and to distinguish the neutrino mass hierarchy. To do this we used muon-like atmospheric 
neutrino events in the GeV range and exploited the matter effects which affect their propagation 
through the Earth. The goal of future multi-Mton detectors must be to lower the energy threshold 
below 10 GeV so that resonant effects can be accessed which would open the road to address many 
important physics questions. 
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A Event numbers 

In this appendix we provide some of the event numbers obtained in our analysis, the details of 
which are given in Sec. [3] in the main text. Tab. [4] shows the event numbers relevant for the 
DeepCore detector set-up and Tab. [5] those for the PINGU-0 configuration, both for 10 years of 
data taking. 



COS 9 


PREM profile 


Ap = +10% 


Ap = -10% 


[-1.0, -0.9] 


2491 (2561) 


2493 (2562) 


2484 (2558) 


[-0.9, -0.8] 


2087 (2261) 


2124 (2273) 


2067 (2253) 


[-0.8, -0.7] 


1363 (1589) 


1419 (1606) 


1315 (1572) 


[-0.7, -0.6] 


665 (813) 


669 (814) 


675 (816) 


[-0.6, -0.5] 


243 (266) 


228 (261) 


261 (273) 


[-0.5, -0.4] 


279 (225) 


274 (224) 


281 (225) 


[-0.4, -0.3] 


919 (857) 


925 (860) 


912 (853) 


[-0.3, -0.2] 


2170 (2133) 


2175 (2136) 


2165 (2130) 


[-0.2, -0.1] 


3814 (3801) 


3816 (3802) 


3813 (3801) 


[-0.1, 0.0] 


5294 (5292) 


5294 (5292) 


5294 (5292) 



Table 4: Total number of muon-like contained events after 10 years in the DeepCore detector 

in the neutrino energy range 10-15 GeV. We show the number of events for each angular bin for the PREM 
density profile (second column) and the cases when the overall normalization is changed by ±10% (third and fourth 
columns). The numbers without (with) parentheses correspond to NH (Iff). For the oscillation parameters we use 
the true values given in Tab. [2] in the main text. 
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cos 9 


PREM profile 


Ap = +10% 


Ap = -10% 


5 - 10 GeV 


10 - 15 GeV 


5 - 10 GeV 


10 - 15 GeV 


5 - 10 GeV 


10 - 15 GeV 


[-1.0, -0.9] 


4031 (4230) 


3649 (3751) 


4114 (4274) 


3652 (3753) 


4164 (4258) 


3638 (3746) 


[-0.9, -0.8] 


4511 (4598) 


3057 (3312) 


4549 (4588) 


3111 (3329) 


4425 (4563) 


3028 (3300) 


[-0.8, -0.7] 


5066 (4815) 


1996 (2327) 


5267 (4856) 


2078 (2353) 


4766 (4725) 


1927 (2303) 


[-U. 1 , -U.Oj 


4832 (5227) 


974 (1191) 


5162 (5379) 


980 (1192) 


4600 (5095) 


989 (1195) 


[-0.6, -0.5] 


5787 (6473) 


356 (390) 


5794 (6518) 


335 (382) 


5869 (6464) 


382 (399) 


[-0.5, -0.4] 


5837 (6196) 


408 (329) 


5689 (6149) 


402 (328) 


6001 (6257) 


411 (330) 


[-0.4, -0.3] 


3260 (3193) 


1346 (1255) 


3193 (3161) 


1355 (1259) 


3320 (3225) 


1336 (1249) 


[-0.3, -0.2] 


1804 (1589) 


3179 (3124) 


1812 (1588) 


3186 (3128) 


1794 (1590) 


3171 (3121) 


[-0.2, -0.1] 


6319 (6202) 


5587 (5568) 


6327 (6203) 


5589 (5569) 


6311 (6201) 


5585 (5567) 


[-0.1, 0.0] 


13513 (13496) 


7754 (7751) 


13514 (13496) 


7754 (7751) 


13513 (13496) 


7754 (7751) 



Table 5: Total number of muon-Hke contained events after 10 years in the PINGU-0 configuration, 

i.e., with two energy bins, 5-10 GeV and 10-15 GeV. We show the number of events for each angular bin for the 
PREM density profile (second and third columns) and the cases when the overall normalization is changed by ±10% 
(fourth to seventh columns). The numbers without (with) parentheses correspond to NH (IH). For the oscillation 
parameters we use the true values given in Tab. [2| in the main text. 



B Note on uncorrelated and correlated errors 

Let us qualitatively explain the effects of uncorrelated and correlated errors on the results. Let 
us consider a total number of bins equal to n, j\T data as the number of events detected in the i-th 
bin and iV* h (A) as the number of predicted events for the parameter set A, with crf tat = \J j\T data 
and <Tj Sys = <7g Nf h (X) being the statistical and systematic error in the i-th bin, respectively, which 
satisfy of at < of 8 . 

For illustration, let us now consider that only the i = 1, ... , m-th bins contribute to the x 2 an d 
that their contribution is approximately the same, i.e., jV data — N- h (X) ~ AiV(A), crf at ~ <x stat and 
<jj Sys ~ <r sys , for i < m. 

In the case of uncorrelated systematic errors, the x 2 (A, c§) function is 



„ [jV data -iV7 h (A)j 2 / f Bten2 \ 

x 2 C\ <*£) = £ VfT + WP * Uf'V+Wv J * 2(A ' 0) ' (33) 

Thus, regardless of how many bins contribute to the \ 2 i the sensitivity gets reduced by a similar 
amount with respect to the case of no systematic errors (cr^ = 0). 



On the other hand, in the case of correlated systematic errors, the X \\ a i) function is 



154- 
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X 2 (A, 0f) = min c 



1 

C7c 



(34) 



where min^ |x 2 (A, cr^, £) j indicates the minimum of x 2 (A, o"£, with respect to £. This is a 
quadratic polynomial in £, so the minimum can be calculated analytically. The result is 



X 2 (A, o*) = x 2 (A, 0) 



EILi #f ^" ata - ^tf)) 



/ s y s 

i±v» ( 



(35) 



Within the approximations we are considering, Eq. 35 can be written as 



x 2 (A,a € ) ~ X 2 (A,0) 



X 2 (A, 0) 



cr 



sys\2 



stat "\ 2 



+ n (a 8 ? 8 ) 2 



cr 



sys\2 



(cr stat ) 2 + n (a s y s ) s 



,i=l 
2/ 



jydata _ jV^(A) 



(7 



stat 



m / ATdata 



(a stat ) 2 + (n - m) (cx sys ) 2 2 - 

( a stat)2 + " ( a s y s)2 X (A, 0) 



r(A, 0) + ^ 

m 



7V* h (A)) (Ar/ ata - A^ th (A)) 



stat ^2 



n 



X 2 (A, 0). 



(36) 



Hence, if most of the bins contribute (m < n), then x 2 (A, o^) <C X 2 (A, 0) and the sensitivity gets 
reduced with respect to the case of no systematic error. However, if only a few bins contribute 
(m <C n), then x 2 (A, 0"^) ~ X 2 (A, 0), so the sensitivity is only slightly affected by the correlated 
systematic errors. 
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